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The first measurements of AxF^ are higher than current theoretical predictions. We investigate the sensitivity of these 
theoretical predictions upon a variety of factors including: renormalization scheme and scale, quark mass effects, higher twist, 
isospin violation, and PDF uncertainties. 



1. Introduction 



Deep inelastic lepton-nucleon scattering experiments 
have provided precision information about the quark 
distributions in the nucleon. However, there has been 
a long-standing discrepancy between the F2 structure 
functions extracted from neutrino and muon experi- 
ments in the small x range. Recently, a new analysis of 
differential cross sections and structure functions from 
CCFR v^-¥e and v^-Fe data was presented; in this 
study, the neutrino-muon difference is resolved by ex- 
tracting the structure functions in a physics model 
independent way.|]lj 

In previous analyses of data,Q structure func- 
tions were extracted by applying a slow rescaling cor- 
rection to correct for the charm mass suppression in 
the final state. In addition, the AxF 3 term (used as 
input in the extraction) was calculated from a lead- 
ing order charm production model. These resulted in 
physics model dependent (PMD) structure functions. 
In the new analysis, jjj slow rescaling corrections are 
not applied, and AxF 3 and F2 were extracted from 
two parameter fits to the data. 

The extracted physics model independent (PMI) val- 
ues for F% are then compared with F% within the 
framework of NLO models for massive charm produc- 
tion; these are found to be in agreement, thus resolving 
the long-standing discrepancy between the two sets of 
data. However, the first measurements of AxF 3 are 
higher than current theoretical predictions. The ob- 
jective of this paper is to investigate the sensitivity of 
AxF 3 upon a variety of factors including renormaliza- 
tion scheme and scale, quark mass effects, higher twist, 
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Figure 1. Comparison of AxF 3 data from CCFR with 
various theoretical predictions. 



isospin violation, and PDF uncertainties. 

In Fig. [l] we have plotted the quantity AxF 3 for an 
isoscalar target computed to order a\. For compari- 
son, we also display data from the CCFR analysis. M 
While there is much freedom in the theoretical calcula- 
tion, the difference between these calculations and the 
data at low Q values warrants further investigation. 
We will discuss and compare the different theoretical 
calculations, and examine the inherent uncertainty in 
each with respect to different input parameters. We 
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will also examine the experimental input, and assess 
uncertainties in this data. 



would still leave a discrepancy between the lowest Q 2 
data points and the theoretical predictions. 



1.1. Measurement of AxF 3 

The structure functions are defined in terms of the 
neutrino-nucleon cross section via: 

d 2 a u ' p 



dx dy 



G F ME W 



7T 



F 2 (l-y)+xF iy 2 ±xF 3 y(l-%) 



(1) 



where Gp is the Fermi weak coupling constant, M is 
the nucleon mass, E u is the incident energy, y = Eh/E u 
is the fractional energy transfer, and Eh is the final 
state hadronic energy. 

The sum of and differential cross sections for 
charged current interactions on an isoscalar target is 
then: 
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where e ~ 2(1 — y)/[l + (1 — y) 2 } is the polarization of 
the virtual W boson. In the above equation, we have 
used the relation: 



2xF 1 (x,Q 2 ) = F 2 (x,Q 2 



l + 4M 2 x 2 /Q< 
l + R(x,Q 2 ) 



(3) 



where R = (Tl/ctt is the ratio of the cross-sections of 
longitudinally- to transversely-polarized VU-bosons, Q 2 
is the square of the four-momentum transfer to the nu- 
cleon, and x = Q 2 /2MEh is the Bjorken scaling vari- 
able. For x < 0.1, R in neutrino scattering is expected 
to be somewhat larger than R for muon scattering be- 
cause of the production of massive charm quarks in the 
final state for the charged current neutrino production. 

Using Eq. (g) and Eq. (§), F 2 and AxF 3 can be ex- 
tracted separately. Because of the positive correlation 
between R and AxF 3 , the extracted values of F 2 are 
rather insensitive to the input R. If a large input R is 
used, a larger value of xF 3 is extracted from the y dis- 
tribution, thus yielding the same value of F 2 . In con- 
trast, the extracted values of AxF 3 are sensitive to the 
assumed value of R, which yields a larger systematic 
error, shown on the data. For R the QCD-inspired fit 
in H is used, but corrected for charged current neu- 
trino scattering using a leading order slow rescaling 
model. This gives precisely the same type of correction 
as a full NLO calculation including the massive charm 
quark, as shown in Ref. 0], but leads to a somewhat 
higher normalization than the perturbative correction. 
It is arguable which prescription for R leads to the bet- 
ter fit to existing data, but the difference between the 
two leads to an uncertainty in AxF 3 of the order of 
the systematic error shown. Clearly a further reduc- 
tion in the assumed value of R (even down to zero), as 
suggested by models including ln(l/x) resummations, 



1.2. Quark Parton Model Relations 

Now that we have outlined the experimental method 
used in the extraction of the F structure functions, it 
is instructive to recall the simple leading-order corre- 
spondence between the F's and the PDF'sfj] 



xF 



x {u + u + d + d + 2s + 2c} 

x {u - u + d - d± 2s =F 2c} (4) 



Therefore, the combination AxF 3 yields: 

AxF 3 = xFg N - xF$ N ~ Ax{s - c} . (5) 

Note that since this quantity involves the parity vio- 
lating structure function F 3l this measurement has no 
analogue in the neutral current photon-exchange pro- 
cess. Also note that since, at leading-order, AxF 3 is 
directly sensitive to the strange and charm distribu- 
tions, this observable can be used to probe the heavy 
quark PDF's, and to understand heavy quark (charm) 
production. We discuss these possibilities further in 
the following subsection. 

1.3. Implications for PDF's 

We have illustrated in Eq. (|J) how AxF 3 is closely 
tied to the heavy quark PDF's. The questions is: given 
the present knowledge base, should we use AxF 3 to 
determine the heavy quark PDF'S, or vice verse. To 
answer this question, we briefly review present mea- 
surements of heavy quark PDF's, and assess their un- 
certainty. 

1.3.1. Tevatron W + Q Production 

The precise measurement of W plus heavy quark 
(W + Q) events provides important information on 
heavy quark PDF's; additionally, such signals are a 
background for Higgs and squark searches. ||,||] 

Unfortunately, a primary uncertainty for W + Q pro- 
duction comes from the heavy quark PDF's. Given 
that AxF 3 is sensitive to these heavy quark PDF's, we 
see at least two scenarios. One possibility is that new 
analysis of present data will resolve this situation prior 
to Run II, and provide precise distributions as an input 
to the Tevatron data analysis. If the situation remains 
unresolved, then new data from Run II may help to 
finally solve this puzzle. In the far future, a neutrino 
experiment from the neutrino factory at a linear col- 
lider would be an ideal tool to measure any neutrino 
structure function. M 



1 To exhibit the basic structure, the above is taken in the limit 
of 4 quarks, a symmetric sea, and a vanishing Cabibbo angle. 
Of course, the actual analysis takes into account the full struc- 
ture.§ 
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Figure 2. a) We display several strange sea PDFs in 
the relevant x range at Q 2 = 5 GeV. b) Same as above 
for the charm sea. 

1.3.2. DIS Di-muon production 

The strange distribution is directly measured by 
dimuon production in neutrino-nucleon scattering.^] 
The basic sub-process is vs — > fi~cX with a subse- 
quent charm decay c — > /i + X'. 

DIS dimuon measurements have safely established 
the breaking of the ST/ (3) flavor symmetry, 



J dx x s{x) 



J dx x [u(x) + d(x)]/2 



< 



(G) 



? 2 =O(10GcV 2 ) 



in the nucleon sea. Still, there remain large uncertain- 
ties for the s-quark distribution in the kinematic re- 
gions relevant for AxF^, even though the (x, Q)-range 
of the CCFR dimuon measurements [0,|l3| is compara- 



ble. CCFR recorded 5044 ^ and 1062 



fj^fj.^ events 



with P w > 9 GeV, > 5 GeV, 30 > E v > 600 GeV, 
and (Q 2 ) = 22.2 GeV 2 . The more recent NuTeV ex- 
periment recorded a similar sample of events, and these 
are presently being analysed through a MC simulation 
based on NLO quark- and gluon-initiated corrections 
at differential levelG3,|l^,f^|. A complete NLO analy- 
sis of this data together with a global analysis may help 
to further constraint the strange quark distribution. 

At present, PDF sets take strangeness suppresssion 
into account by imposing the constraint in Eq. (|^) as 
s(x) — k[u(x) + d(x)}/2 at the PDF input scale Qo ~ 
1 GeV, p6|,p^| or by evolving s(x) from a vanishing 
input at a lower scale. flcl The residual uncertainty can 
be large, as can be seen from the collection of strange 
seas in Fig. ^. 

In the following, we will take the uncertainty in the 
strange PDF into account by relaxing the experimen- 
tal constraint in Eq. (](]) as implemented in the MRST 
partons.|l(| This will be discussed in Sec. |3.2|. 



2 Presently, there are a number of LO analyses, and one NLO 
ACOT analysis. Hfl. Results of a recent LO analysis by NO- 



1.3.3. Charged and Neutral Current DIS 

The strange distribution can also be extracted indi- 
rectly using a combination of charged (F%) and neutral 
(FJf) current structure functions; however, the system- 
atic uncertainties involved in this procedure make an 
accurate determination difficult. |l| The basic idea is to 
use the (leading-order) relation 

(c + c) + 

18 " ' 



(7) 



2 ±<J L 5 9 + 9 

to extract the strange distribution. Here, q + q rep- 
resents a sum over all quark flavors. This method is 
complicated by a number of issues including the xF% 
component which can play a crucial role in the small- 
x region — precisely the region where we observe the 
discrepancy. From the corresponding relation: 



:Fi 



r 2 



1 



18- --12^ (8) 
we see that these problems are not independent; how- 
ever, this information, together with the exclusive 
dimuon events, may provide a more precise determi- 
nation of the strange quark sea, and help to resolve 
our puzzle. 

Prior to the DIS dimuon data, the 1992 CTEQ1 
analysis [^9| found that a combination of NC struc- 
ture functions from NMCp0[ and the physics-model- 
dependent charged current structure functions from 
CCFR H seem compatible with approximate SU(3)- 
symmetry, i.e., k ~ 1 in Eq. (^|). Recent dimuon mea- 
surements now exclude an St/(3)-symmetry k = 1. 



We shall explore the effect of k on AxF% in Sec. [3.2] 



MAD [12[ are inline with these experiments. 



2. Dependence of AxF$ on input parameters 

We now systematically investigate the sensitivity of 
the theoretical predictions of AxF^ upon a variety of 
factors including: renormalization scheme and scale, 
quark mass effects, higher twist, isospin violation, and 
PDF uncertainties. 

To simplify this analysis, we first examine the influ- 
ence of these factors on the LO expression: Aici 7 -} = 
4x{s(£, fi) — c(x, /j,)}; after using this as a "toy model," 
we will then return to the full NLO calculation in 
the next section.^] For most variables, the simpli- 
fied LO is sufficient to display the general behavior 
of the full NLO result. There are two exceptions: 1) 
the scheme dependence, and 2) the PDF dependence. 
These factors depend on the interplay of both the 
quark-initiated LO contributions as well as the NLO 
gluon-initiated contributions. For this reason, we will 
postpone discussion of these effects until the following 
section. 

2.1. Charm Mass 

We start by examining the effect of the charm mass, 
m c , on AxF^. In Fig. [j[ we plot the leading-order 

3 Note that a possible s ^ S asymmetry JH| would average out 
in AxF-j. 



3 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



4 


» 


m c =0.0 












X=0.015 



x=0.015 



5 10 20 

Q 2 (GeV 2 ) 




5 10 20 

Q 2 (GeV 2 ) 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



x=0.08 








' m c =0.0 


r 2 


S 10 


m c =1.8 

20 50 1 00 



Q 2 (GeV 2 ) 

Figure 3. Variation of the leading-order expression 
AxF% — 4x{s(£, /i) — c(x, /i)} on the charm mass, m c , 
plotted vs. fi. We take values of m c in the range 
m c = [0,1.8] GeV using steps of 0.2 GeV for three 
choices of x. We use £ = x(l + m 2 /Q 2 ), and li 2 = 
(Q 2 + 1.6 2 GeV 2 ). 



Q 2 for 



expression AxF^ = 4x{s(^,^) 
values of m c in the range = [0, 1.8] GeV using steps 
of 0.2 GeV. Here, we define £ = x(l + rn 2 /Q 2 ) which 
is a "slow-rescaling" type of correction [g^,|23|] which 
(crudely) includes mass effects by shifting the x vari- 
able. "Slow rescaling" naturally arises at LO from 
the single-charm threshold condition(VF > m c ) and 
is required at NLO for consistent mass factorization 

nil- 

Note, the result of this correction is most significant 
at low Q 2 . To isolate the kinematic m c dependence, we 
have chosen the scale li 2 — Q 2 + 1.6 2 GeV 2 . (We will 
separately vary the scale in the following subsection.) 

Note that this exercise is only altering the charm 
mass in one aspect of the calculation; to be entirely 
consistent it would be necessary to obtain parton dis- 
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Figure 4. Variation of the leading-order AxF^ on 
the renormalization scale, li. We have chosen /x 2 = 
{Q 2 ,Q 2 + m 2 ,P 2 }. 



tributions (particularly the charm quark) using fits 
with different charm masses. For charm masses in the 
range 1.2 to 1.8 GeV, this will be a small effect; for 
charm masses below ~ 1.2 GeV, (which is below the ex- 
perimentally allowed range |25|, pq ) , such issues become 
important and the curves of Fig. [| will be modified, 
i.e. lower m c will lead to longer evolution for charm, a 
larger charm distribution, and a lowering of the curves 
in Fig. | 

While taking m c — > does raise the theoretical 
curves in the regions where we observe the discrepancy 
(namely, the low Q region), varying m c even within the 
wide range of [1.2, 1.8] GeV (lower 4 curves) does not 
give us sufficient flexibility to match either the shape 
or normalization of the data. 

2.2. ix Scale 

Next, we investigate the variation 
of AxF^ = Ax{s(^,fj.) — c(x,[i)} with the renormal- 
ization/factorization scale ll. We use three choices of 
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Q (GeV) 

Figure 5. a) We display P™ ax and Q vs. x for 
E=200 GeV and y = {0.2,0.4,0.6,0.8}. b) Pp ax vs. 
Q for x = {0.3,0.1,0.03,0.01}. 

the p? scale: 

• M 2 = Q 2 , 

• /j 2 = Q 2 + m 2 , 



p2 



Of course, AxF 3 with /i 2 scales of Q 2 and Q 2 + m 2 
differ only at lower values of Q 2 ; AxF% with /j 2 scale 
of PT ma7 . is comparable to Q 2 and Q 2 + m 2 at larger a;, 
but lies below for smaller x. The scale choice Q 2 + m 2 
leads to an improvment over Q 2 by providing a lower 
bound on /x to keep the scale in the perturbative region. 

The choice of Pr meu , is motivated, in part, by some 
observations by Collins. [^7| To display the relationship 
between Pr max and Q, in Fig. || we plot both Pr max 
and Q vs. x for 4 choices of y; note]l3]| that the x- 
dependence of Pr max is opposite that of Q. While tying 
the scale choice of \i to PT max has some interesting 
intuitive interpretations,^ for small x this scale clearly 
becomes too large for the relevant physics. In any case, 
it cannot help us with our AxF3 problem as the scale 
choice of Pr ma:c moves the theory curves away from the 
data. 

2.3. Higher Twist 

We now illustrate the potential effects due to higher 
twist contributions. We parameterize such contribu- 
tions by multiplying the leading-twist terms by a cor- 



4 Technically, the interpretation is in terms of the characteristic 
Pt of the partonic subprocess, but this is unobservable; therefore 
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Figure 6. Variation of the LO AxFs, on higher twist. 
The functional form is (1 + D^{x) / Q 2 ) where we take 
D 3 (x) in the range D 3 (x) = [0,10] GeV 2 in steps of 
1 GeV 2 . 

rection factor that will increase at low Q 2 . More specif- 
ically we use:[f28[ 



F i HT (x,c?)=Fr(x ) cr)[i+ 



LT i 



Di{x) 

Q 2 



(9) 



is used instead. 



where i = {1,2,3}. We vary D 3 (x) over the range 
[0,10] GeV 2 in steps of 1 GeV 2 . (For the purposes of 
our simple illustration, it is sufficient to take D 3 (x) to 
be independent oi x.) We will find that this range is 
well beyond what is allowed by experiment; however, 
we display this exaggerated range to make the effect of 
the higher twist contributions evident. 

To normalize this choice with the allowable range 
consistent with data, we compare with the MRST 
higher twist fit which extracted a limit on the function 
£>2 (x) . We see from the table of Z?2 (x) that the allowed 
contribution from the higher twist terms is quite small. 
In addition, we note that the sign for D2(x) obtained 
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for the relevant small- x region tends to be negative- 
exactly the opposite sign that is needed to move the 
theory toward the data. 

As we have shown constraints on D2 (x) , the obvious 
question is should we expect D3 (x) to be substantially 
different? A calculation of the power corrections using 
renormalons p9[ suggests that while both D\(x) and 
D^(x) are of the same order of magnitude as D2{x) 
and have similar x-dependence, that D\{x) and D^{x) 
are even more negative than D2(x); again, this trend 
would move the theory farther from the data. There- 
fore, we conclude that that D2(x) represents a conser- 
vative limit for the D^{x) power corrections. 

Examining Fig. ^, we note that it would take an 
enormous higher twist contribution to bring the nor- 
malization of the curves in the range of the data points, 
and still the shape of the Q 2 -dependence is not well 
matched; hence, we conclude that this is not a com- 
pelling solution. 
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Table 1 



Values of the higher-twist coefficient D2(x) extracted 
from the LO, NLO and NNLO fits. Units are GeV 2 . 
Table taken from Ref. |FJ. 



2.4. Isospin Violations 

The naive parton model identity in Eq. (||) is mod- 
ified if the full (non-diagonal) CKM structure, NLO 
QCD radiative corrections, and QCD based charm pro- 
duction are taken into account. This expression is 
also modified even in leading- order if we have a vio- 
lation of exact p <-> n isospin-symmetry, (or charge 
symmetry); e.g., u n (x) ^ d p (x). In deriving Eq. (||), 
isospin-symmetry was necessary to guarantee that the 
u, d-contributions cancel out in the difference, thereby 
leaving only the s(x) and c{x) contributions.^] 

5 Note we have not investigated shadowing corrections. The 
CCFR data on Fe is converted with an isoscalar correction, and 
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Figure 7. Impact of a one-parametric charge symme- 
try violation (CSV) toy model of Eq. © and Eq. (0) 
on the LO AxF 3 of Eq. ©. 6 is varied over the range 
-6 = [0, tt/2] in steps of (tt/20). 

The validity of exact charge symmetry (CS) 
has recently been reexamined. |]3l| Residual u,d- 
contributions to AxF 3 from charge symmetry violation 
(CSV) would be amplified due to enhanced valence 
components {u v (x), d v (x)}, and because the d — > u 
transitions are not subject to slow-rescaling correc- 
tions which strongly suppress the s — > c contribution 
to AxF%, as illustrated in Fig. |^. 

We will examine possible contributions to AxF% by 
considering a "toy" model to parameterize the CS- 
violations. In isospin space, we can parameterize a 
general transformation as a rotation: 

k)^ sv = N q J2R qq >(e)W) P (10) 

the corresponding uncertainties are included in the data. O 
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where R is a rotation matrix, and N q is the normal- 
ization factor. For example, in this model the it- 
distribution in the neutron would be related to the 
proton distributions via the relation: 



,csv 



N 2 [cos 2 (9)u p (x, Q 2 ) + sin 2 (0)d p (x, Q 2 )] (11) 

For 9 = 7r/2, we recover the symmetric limit 
u p (x,Q 2 ) = d n (x,Q 2 ). Note that we define the nor- 
malization N u such that this model preserves the sum 
rule: 



-CSV 



(x,Q 2 ) = l 



(12) 



Note that Eq. (|10|) should not be considered a serious 
theory,^] but rather a simple one-parameter (9) model 
which is flexible enough to illustrate a range of CSV 
effects for AxF 3 . 

In Fig. ^ we vary 9 over the its maximum range 
[0,7r/2] in steps of n/20. The exact charge symmetry 
limit (9q — 7r/2) of Eq. ( pd| ) corresponds to the lowest 
curve in Fig. [?]. Note the effect of the CSV contribution 
monotonically increases as 9 deviates from the charge 
symmetry (CS) limit 9q = tt/2. From the plot we 
observe that a violation of 9o — 9 > (3/10)(tt/2) is 
required to bring the theory in the neighborhood of 
the data. 

At the relevant x values of interest for AxF 3 , this 
translates (via Eq. ( [TT| ) and its analogue for d n ) into 
a ~ 20% symmetry violation for u n , and a ~ 10% 
symmetry violation for d n . Specifically, 

d p (x) 



u^{x) 

u p {x) 
d^(x) 



1.2 



0.9 



(13) 
(14) 



Since the electroweak couplings are flavor-independent, 
( V q — A q — 1 V q) , AxF 3 is in principle insensitive to a 
re-shuffling of the CSV contributions between u n and 
d n . In particular, if we define the shift due to CSV in 
AxF 3 as 5AxF 3 , we find: 



SAxF 3 ee - (Sd + 5d - Su - Su) 
where 



(15) 



(16) 
5|) is eval- 



Su{x,Q 2 ) = u^{x,Q 2 )-d p (x,Q 2 ) . 

The expression for SAxF 3 (x, Q 2 ) of Eq 
uated at Q 2 = 4 GeV 2 from Eq. fl) with 9 a - 9 = 
(3/10)(tt/2) and is plotted in Fig. §|] 

Although this level of isospin violation certainly im- 
proves the description of AxF 3 , it is necessary to con- 
sider precisely what level of violation is actually al- 
lowed by other experimental data. For instance, it 

6 E.g., it does strictly speaking not commute with evolution. 
7 Again, the detailed x-shape from Eq. (|l]]) should not be taken 
as a serious model prediction. 
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in Eq. @ and Eq. (Ill) for 9 - 9 = 3(tt/20). 



has previously been suggested that the discrepancy be- 
tween F 2 from neutrino and muon data itself may be 
due to isospin violation. [pT| This type of violation re- 
quired that u be very different to d in the region of 
interest. 

The measurement of the lepton charge asymmetry in 
W decays from the Tevatron |3^,|33| places tight con- 
straints on the up and down quark distributions in the 
range 0.007 < x < 0.24, constraining them to be ap- 
proximately as specified in the parton sets obtained by 
the global analyses. While only strictly telling us about 
parton distributions in the proton, this data rules out 
an isospin violation of this type to about 5% as demon- 
strated in 

However, there are other strong constraints on 
isospin violation. For example, we note that while the 
toy model above leaves the neutron singlet combina- 
tion q + q invariant at the < 2% level in the region 
x e [0.01; 0.1], it would lower the NC observable: 



1 



(u + u) + -(d + d) 



(17) 

a;e[0.01;0.1] 

by about 10%. An effect of this size would definitely be 
visible in NMC F^ / F% data which has an uncertainty 
of order a few percent in this kinematic region, and 
acts as a major constraint, [pof 

At this point, one could play clever games to evade 
the constraints of specific experiments. For example, a 
re-shuffling of CSV contributions between the individ- 
ual Sq, 5q in Eq. ([If]) according to: 



xdn 



xq n 



- -6AxF 3 
5 

+ ^5AxF 3 



(18) 



analogous 



would keep Eq. (17) invariant. However, this would in 
turn raise the CC observable Eq. (|l5| ) by > 5%; though 
it would help to explain the excess in 5xF 3 , it would 
spoil the new-found compatibility between neutral cur- 
rent and charged current data. 
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Figure 9. LO and NLO calculations for AxF 3 vs. Q 2 
in 3 and 4 flavor schemes for three values of x. The 
dotted curve (f=3, a°) corresponds to the 3- flavor LO 
result. The dashed curve (P=4, a®) corresponds to the 
4-flavor LO result. The solid and dot-dash curves cor- 
respond to the 3 and 4 flavor NLO QCD calculations, 
respectively. 



In addition, there are also fixed-target Drell-Yan ex- 
periments |54|,|3^[ such as NA51 and E866 which pre- 
cisely measure d/u in the range ]3^| 0.04 < x < 0.27, 
and are also sensitive to isospin violating effects. 

We therefore conclude that the many precise data 
sets which constrain different combinations of the 
PDF's probably leave no room for CSV contributions 
of the magnitude necessary to fully align the theory 
curves with the AxF 3 data. However, we should add 
the caveat that an exhaustive investigation of the in- 
terplay of these different data sets and their influence 
on AxF 3 will only be possible within a global PDF 
analysis 

3. NLO calculation of AxF 3 /2 

Now that we have used the leading-order expression 
for AxF 3 to systematically investigate dependence of 
this observable on various parameters, in this section 
we now turn to the full NLO calculation. 

3.1. Contributions to the NLO Calculation 

In Fig. [| we have plotted the LO and NLO calcula- 
tions for AxF 3 vs. Q 2 on an isoscalar target in 3 and 
4 flavor schemes. The 3-flavor LO calculation (f=3, 



a°) involves primarily the strange quark contribution, 
s(x), as the charm distribution is excluded in this case. 
When the higher order terms are included (f=3, a\), 
this result moves (substantially) toward the predictions 
of the 4-flavor scheme. We note that while the 3-flavor 
LO calculation (f=3, oP s ) appears consistent with the 
data, we cannot take this result as a precise theoretical 
prediction as this simplistic result is highly dependent 
on scheme and scale choices; a result that is verified by 
the large shift in going from LO to NLO. 

The pair of curves in the 4-flavor scheme (using the 
CTEQ4HQ distributions) nicely illustrates how the 
charm distribution c(x,/j, 2 ) evolves as \n(Q 2 /m 2 ) for 
increasing Q 2 ; note, c(x,/j, 2 ) enters with a negative 
sign so that the 4-flavor result is below the 3-flavor 
curve. For the scale choice, we take /i = \JQ 2 + m 2 . 
While the scale choice /z = Q is useful for instruc- 
tive purposes such as demonstrating the matching of 
the 3- and 4-flavor calculations at /i = Q = m c , the 
choice fj, = \/ Q 2 + m 2 is more practical as it provides 
a lower bound on /i which is important for the PDF's 
and a a {n). (C/., Sec. [|J], and Ref. §§.) 

Additionally we note the stability of the 4-flavor 
scheme in contrast to the 3-flavor scheme. The shift 
of the curves when including the NLO contributions is 
quite minimal, particularly when compared with the 3- 
flavor result. p7|,|3^] This suggests that organizing the 
calculation to include the charm quark as a proton con- 
stituent can be advantageous even at relatively low val- 
ues of the energy scale. 

3.2. PDF Uncertainties: s(x), ... 

In Fig. [HJ we show the variation of NLO calculation 
of AxF 3 on the strange-quark PDF. To obtain a real- 
istic assessment of the s(x) dependence, we have use 
the NLO calculation with PDF's based on the MRST 
set which are re- fit with the value of k = 2s/ (u + d) 
constrained to be k = {0.50,0.78,1.00}. Note, by re- 
fitting the PDF's with the chosen value of k we are as- 
sured to have an internally consistent set of PDF's with 
appropriate matching between the quarks and gluon, 
and with the sum-rules satisfied^] 

The choice k = 0.50 is in line with the many experi- 
mental determinations of n, c/., Table. |^; as expected, 
this prediction lies farthest from the data points. 

The choice k — 0.78 is taken as an extreme upper 
limit given the experimental constraints; actually, in 
light of the results of Table. |[ this is arguably be- 
yond present experimental bounds. This prediction is 
marginally consistent at the outer reach of the system- 
atic + statistical error bars. 

Finally, we take an SU(3) symmetric set (k = 1) 



8 Note, re is certainly Q-dependent, and the values for re quoted 
above correspond to the Qq of the evolution. While re compares 
the integral of s(x) to the seq-quarks, there is also the possibility 
of an x-dependent variation. This has been studied in the fits of 
the strange-sea J9| [ltj, [^J ; we shall find that such subtle effects 
can play no role in resolving the AxF'j issue. 
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Figure 10. Variation of Axi^ on the strange-quark 
PDF. The NLO calculation is used with three sets 
of PDF's. These PDF sets are re-fit based on the 
MRST set, and the value of k is constrained to be 
k = {0.50,0.78,1.00}. 

purely for illustrative purposes. It is interesting to note 
that even this extreme value is still below the central 
value of the data points at the higher x values. 

In conclusion we note that increasing the strange 
quark distribution does succeed in moving the theory 
toward the data; however, our consistent NLO analysis 
presented here suggests that the we have only limited 
freedom to increase s(x), and that this alone is not suf- 
ficient to obtain good agreement between theory and 
data. 

3.3. Scheme Choice 

In our final section, we present the best theoreti- 
cal predictions presently available to demonstrate the 
scheme dependence of AxFs(x, Q). Specifically, in 
Fig. [ll] we show predictions for: 

• NLO FFS GRV.Jl| 

• NLO VFS TR calculation. [0 

• NLO VFS ACOT calculation with CTEQ4 

pdf's. ©mm 

• NLO VFS ACOT calculation with CTEQ5 
PDF's. QUE! 

All these calculations use NLO matrix elements, and 
are matched with appropriate global PDF's which are 



Experiment 


Order 


K 


Rcf. 


CDHS 


LO 


47 ± 08 ± 05 




IN 
E| 




FMMF 


LO 


n A i +0.075 +0.103 
— u.u/o — u.uoy 


I 


40 


1 


CHARM II 


LO 


n oq+O.OY +0.07 
U - OJ -0.06 -0.07 




11 


1 


CCFR* 


NLO 


477+U.U46 +0.023 
u -^' '-0.044 -0.024 


I 


13 


1 


CCFRf 


NLO 


/lfi«+ () - lM;i + U.U24 
U.1UO_Q.0 46 -0.025 


I 


13 


1 


CCFRf 


LO 


0.373+^ ± 0.018 
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NOMAD 


LO 


n 40+u.uy +0.17 

u -^°-0.07 -0.12 


1 


12 
10 


L 

1 


NuTeV 


LO 


0.42 ±0.07 ±0.06 


I 



Table 2 

Next-to- leading-order and leading-order fit results. Er- 
rors are statistical and systematic. This table is dis- 
played to estimate the upper limits allowed by experi- 
ment; a full comparison must take into account scheme 
and scale choices, and the shape parameters. 
* Collins- Spiller Fragmentation, 
t Peterson Fragmentation. 



fitted in the proper scheme. 

The first observation we make is how closely these 
four predictions match, especially given the wide vari- 
ation displayed in previous plots such as Fig. ^. In 
hindsight, this result is simply a consequence of the fact 
that while different renormalization schemes can pro- 
duce different results, this difference can only be higher 
order. [] Thus, the difference between these curves is 
indicative of terms of order a 2 s which have yet to be 
calculated.^ When terms of order are included, 
the span of these predictions will be systematically re- 
duced to order a^ +1 . 

In Fig. [H], we note the very close agreement among 
the VFS calculations, particularly the TR calcula- 
tion and the ACOT calculation with CTEQ4 PDF's. 
The ACOT calculation with the two CTEQ curves 
show primarily the effect of the charm distribution, 
as CTEQ4 uses m c = 1.6 and CTEQ5 uses m c = 1.3. 
The GRV calculation shows the effect of using yet a 
different scheme, in this case a FFN scheme, with its 
appropriately matched PDF. Were we to use MRST or 
CTEQ PDF's, the spread of these theory curves would 
decrease; however, this would most likely represent an 
underestimate of the true theoretical uncertainties aris- 
ing from both the hard cross section and PDF's.[^] 

While we consider it a triumph of QCD that dif- 
ferent schemes truly yield comparable results (higher 
order terms aside), we should be cautious and note 
that the spread of these curves can only underestimate 
the true theoretical uncertainty. Note that GRV has a 



9 To be precise, different renormalization schemes can differ by 
i) terms of higher order in the perturbation series, and ii) terms 
of higher twist which do not factorize. |l4|, S 
10 For asymptotic results at order a^, see Rcf. [pTJ 
11 The computation of PDF errors is a compL 
some recent approaches to this topic see: Refs.[ 



subject. For 



9 




Figure 11. Variation of AxF^ on the renormaliza- 
tion scheme. All curves use NLO calculations, and 
appropriately matched PDF's. We note that the theo- 
retical predictions are consistent within the theoretical 
uncertainty — as they should be. 



rather different strange distribution due to a different 
philosophy of obtaining this distribution rather than 
due to a different scheme. 

4. Conclusions and Outlook 

Comprehensive analysis of the neutrino data sets can 
provide incisive tests of the theoretical methods, par- 
ticularly in the low Q 2 regime, and enable precise pre- 
dictions that will facilitate new particle searches by 
constraining the PDF's. This document serves as a 
progress report, and work on these topics will continue 
in the future. 

Theoretical predictions for AxF^ undershoot pre- 
liminary fixed target data at the ~ 1 er-level at low 
x and Q. The neutrino structure function AxF% is 
obviously sensitive jl7) to the strange sea of the nu- 
cleon and the details of deep inelastic charm produc- 
tion. A closer inspection reveals, however, consider- 
able dependence upon factors such as the charm mass, 
factorization scale, higher twists, contributions from 
longitudinal polarization states, nuclear shadow- 
ing, charge symmetry violation, and the PDF's. This 
makes AxF^,(x,Q 2 ) an excellent tool to probe both 
perturbative and non-perturbative QCD. 

We have explored the variation of AxF^,(x,Q 2 ) on 
the above factors and found none of these to be capa- 



ble of resolving the discrepancy between the data and 
theory. 

Although we have not eliminated the possibility that 
the entire set of parameters conspires to align the the- 
ory with the data, we have demonstrated this to be 
extremely unlikely. Of course, a definitive answer can 
only be obtained by a global analysis which combines 
the neutrino data for dimuons, AxF$, F^PMI), and 
e 1 ■ 

As the situation stands now, this AxF 3 (x, Q 2 ) puz- 
zle poses an important challenge to our understanding 
of QCD and the related nuclear processes in an impor- 
tant kinematic region. The resolution of this puzzle 
is important for future data analysis, and the solution 
is sure to be enlightening, and allow us to expand the 
applicable regime of the QCD theory. 
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